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Abstract 

We predict properties of triplet excited states in single-walled carbon nanotubes (CNTs) using a 
time-dependent density- functional theory (TD-DFT). We show that the lowest triplet state energy 
in CNTs to be about 0.2-0.3 eV lower than the lowest singlet states. Like in vr-conjugated polymers, 
the lowest CNT triplets are spatially localized. These states show strong optical absorption at 
about 0.5-0.6 eV to the higher lying delocalized triplet states. These results demonstrate striking 
similarity of the electronic features between CNTs and vr-conjugated polymers and provide explicit 
guidelines for spectroscopic detection of CNT triplet states. 
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Semiconductor or metal-like electronic features of single-walled carbon nanotubes (CNTs) 
result from a delocalized vr-electron system confined by a quasi-one-dimensional structure 
of the material^. Understanding, probing and controlling fundamental electronic proper- 
ties of CNTs are the keys opening many exciting nanotechnology applications as nanoelec- 
tronic and nanophotonic devices^i^. In the past 2-3 years we have observed a revolution 
in experimental studies related to carbon nanotube photophysics. Experimental^'^'^'^ and 
theoretical^i^ii^'iiii^'i^'i^ work has revealed the importance of excitonic affects and electron- 
electron correlations enhanced by one- dimensional confinement conditions. Moreover, re- 
cent research has demonstrated substantial electron-phonon coupling (vibrational effects) in 
CNTs as weiy^ii^>i^ii^>i2,. AH these phenomena are typical features of quasi-one- dimensional 
materials^ such as conjugated organic^ and organometallic^^ polymers, and mix-valence 
chains2^. In particular, recent studies reveal many common electronic features observed in 
spectra of CNTs and conjugated polymers^i2iii^'2ii2^. Figure [1] schematically shows the es- 
sential electronic states contributing to photophysics of these materials providing common 
comparison baseline. We use shorthand labeling of the relevant ground (.So) and excited 
states {Si, Ti, . . .), which is common for molecular physics, and field-specific notations. 

So far, extensive amount of work in nanotube science has been done on the excitonic 
properties of the singlet manifold only^^'^^'^'^^^^. Triplet states in CNTs have not been 
detected experimentally yet. Properties of triplet states in CNTs should affect a number of 
fundamental physics phenomena. Singlet-triplet splitting in low-dimensional materials is a 
measure of electronic correlation strength and exchange effects, similar to that of exciton 
binding energy shows. It is not known how relaxation of photoexcitations occurs in CNTs^ii^. 
Carbon nanotubes are weakly-emissive materials. Recent studies attribute this property to 
the 'dark' singlet excitons below the optically allowed staiesr-^^. It is remains to be seen 
if the triplet states play any substantial role in the non-radiative decay of photoexcitations. 
Note that in spite of weak spin-orbit coupling, a substantial fraction of photoexcitations 
undergoes intersystem crossing to triplets in such related systems as fuUerenes. Likewise, 
photoexcitations to the higher lying CNTs excitons (e.g. E22, E33) generate unbound 
electron-hole pairs, which can recombine back to excitons formally at a statistical ratio 
(three triplet states are formed per one singlet). Subsequently, CNTs triplet states can play 
even more important role for relaxation processes from the higher excited states^^i^i or for 
the dynamics of electron and hole carriers in transport processes^"^. Moreover, triplet states 
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may be involved in a number of other important fundamental phenomena such as impact 
ionization and electro-luminescence. 

In this letter we report extensive computational modeling of excited states of CNTs using 
first principles based time-dependent density functional theory (TD-DFT). First of all, we 
ensure that the hybrid DFT used is quantitatively accurate for both singlet and triplet states 
in conjugated polymers and singlet states in CNTs. This warrants quantitative prediction 
of the lowest triplet states in CNTs to be about 0.2-0.3 eV lower than the lowest singlet 
states, and the triplet state absorption pronounced at about 0.5-0.6 eV. We further compare 
properties of calculated singlet and triplet states in polymers and nanotubes using real- 
space transition density matrix analysis. This establishes similarity in electronic properties 
of these materials and guide design of new experimental probes of triplet states in CNTs. 

We focus on (7,6) semiconductor tube segment and poly(p-phenylene vinylene) (PPV) 
oligomer shown in Fig. [2] (top). Reference PPV compound is one of the most explored 
and understood luminescent polymer—, and (7,6) tubes are common species in experimental 
samples. These molecular systems have a finite length of 10 nm and comprise 2 (tube) and 
16 (polymer) repeat units. The finite size one-dimensional systems, when their lengths are 
significantly larger than the characteristic exciton sizes, are expected to reproduce well the 
properties in the infinite-size limit. Calculations of conjugated oligomers provide a standard 
example of such an approacb^ii^^. Unsaturated chemical bonds at the CNT ends have been 
saturated with hydrogen atoms to remove mid-gap states caused by dangling bonds. The 
Austin Model 1 (AMI) Hamiltonian^^ has been further used to obtain optimal geometries of 
both molecules. This method works very well for geometry optimizations in a broad variety 
of vr-conjugated molecular systems. Geometry optimizations conducted with other methods 
would lead to small and uniform red (DFT geometries) or blue (Hartree-Fock geometries) 
shifts of CNT excitation energies (see supplemental materials). Based on the optimized ge- 
ometries, we calculate the corresponding TD-DFT excited state structures up to 25 lowest 
excited states in singlet and triplet manifolds using the Gaussian 03 package^. TD-DFT 
is currently a mainstream approach for quantitative modeling of optical responses in large 
molecules^-. We use B3LYP and PBEIPBE hybrid functionals, containing 20% and 25% 
of the orbital exchange, respectively. These functionals are subjectively considered to be 
the most accurate in computational chemistry reproducing well electronic excitations in 
many materials^. We emphasize the necessity of using hybrid DFT to account for excitonic 
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effects^, which are important for molecules in question. Furthermore, the delicate interplay 
between energetics of singlet and triplet states and respective exciton binding energies are 
extremely sensitive to this exchange component^. Both minimal valence (ST0-3G) and 
extended (6-31G) basis sets were used for all calculations. Currently such TD-DFT (e.g. 
B3LYP/6-31G) calculations became standard in the molecular modeling. However, com- 
puting many excited states in (7,6) tube with more than 1,000 atoms (about 10,000 basis 
functions) is extremely numerically demanding and memory intensive task. 

Table [T] summarizes our main computational findings. First of all, we observe several 
generic trends expected from the methods used. Minimal basis set (ST0-3G) is clearly 
lacking 'room' for electronic delocalization. Subsequently, its excitation energies are shifted 
to the blue compared to the larger 6-31G basis set results. We note that the difference is 
smaller in the nanotube case. Likely, circumferential dimension of CNT somewhat compen- 
sates the reduced basis set size. Further increase of the basis set (e.g. 6-31G*) will likely 
have a minor effect on the CNT excited state energetics leading to red-shifts of excitation 
energies up to 0.1 eV across the board (see supplemental materials). Compared to B3LYP, 
PBEIPBE singlet excitation energies are blue-shifted. This trend is reversed for the triplet 
states. This phenomenon is expected as well: larger fraction of the orbital exchange results 
in increased shifts of singlet and triplet states up and down the spectrum, respectively^ii^^. 
Overall, an accuracy of the calculated excitation energies improves when going from the left 
to the right in table [l PBE1PBE/6-31G or B3LYP/6-31G values agree well (within 0.1-0.3 
eV) with available experimental data across the board. Note that the PPV oligomer was 
calculated in a planar conformation. Accounting for a torsional disorder typically present 
in experimental samples would shift the calculated singlet frequencies to the blue by ~ 0.2 
eV improving an agreement with experiments^. Mainly the trend between computed and 
experimental values is consistent for all excited state of a given material. 

To analyze the excited state properties we use the correspondent transition density ma- 
trices, representing the electronic transition between the ground state and an electronically 
excited states^'^'^. The relevant excitonic states of all molecules are shown in Fig. O The 
matrix diagonal and off-diagonal sizes characterize the distribution of an excitonic wave- 
function over the chain, namely center-of-mass and distance between particles, respectively. 
The 5*1 transition is shown in the first row. Due to an enhanced excitonic character, the 
lowest exciton 'collects' most of the oscillator strength from the band and, subsequently, 
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transition 5*0 — is strongly optically allowed in both PPV and (7,6) tube. The contour 
plots illustrate that the center of mass of photoexcited electron-hole pair is delocalized over 
the entire chain (diagonal in the plot) representing typical vr — vr excitations. According to 
the color code, the exciton size (maximal distance between electron and hole) measured by 
off-diagonal extent of the non-zero matrix area is about 3 nm in PPV and 4 nm in nanotube 
(50% drop of the wavefunction) and about 5 nm in PPV and 7 nm in nanotube (90% drop 
of the wavefunction). Si state is the lowest singlet in the PPV oligomers, which exemplifies 
a typical case for all photoluminescent polymers. The situation is different in CNTs: our 
calculations result in several optically forbidden (or nearly forbidden) exciton states denoted 
as Dll in Fig. [H slightly below (up to 0.1 eV) the allowed Si state. This relative ordering 
of 'bright' and 'dark' states may be responsible for the poor fluorescence efficiency of the 
CNTs^. These states (not shown) have approximately the same delocalization properties as 
Si state displayed in Fig. [2] (top). 

Sn (or mAg) state is the next essential state in the singlet excitonic manifolds of both 
polymers and nanotubes. This state shows up as a major peak in both photoinduced (PAl 
band) and two-photon absorptions of nanotubes as evidenced by two recent experimental 
studies^'Si. The group theory based on the band model does not predict the selection 
rules needed to explain the two photon experiments in CNTs^^. Yet, in the finite molecule 
calculations presented here, Sn state characteristically appears as a state with a strong 
transition dipole moment from 5*1 state. It is, however, strictly forbidden in the ground- 
state absorption. Symmetry notations lAg, IBu, mAg for 5*0, Si, and Sn states, respectively, 
common in the polymer's photophysics, rationalize analogous appearance of Sn state in 
nanotubes as well (a detailed joint experimental/theoretical study appeared recently^^). 
In the real-space analysis, Sn state is much more delocalized compared to 5*1 state with 
significant spatial separation between an electron and a hole (second row in Fig. [2]), thus 
representing weakly bound exciton. Subsequently, the energetic separation between Si and 
Sn state has been used as a lower bound estimate of the 5*1 excitonic binding energy in 
photoluminescent polymers^^ (about 0.6-0.8 eV) and, recently, in CNTs^"^ (about 0.3-0.4 
eV). Above Sn there are several excited states known in polymer's photophysics as nB^ and 
kAg, which represent non-interacting electron- hole pairs at the continuum of the excitonic 
band. Similar states have been recently observed in CNTs spectra as well^i. AH singlet states 
discussed above constitute an elegant essential state model applicable for both polymers and 
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nanotubes . 

Spectroscopic study of triplet states in materials with small spin-orbit coupling is a chal- 
lenge for experiment. Yet, the triplet energies has been measured in a broad range of different 
7r-conjugated polymers^I. This task remains a problem for CNTs. Comparing the relative 
energies of the singlet and triplet gaps provides an alternative estimate of the strength of 
electron-electron correlations, and, particularly, electronic exchange effects. Typically in 
one-dimensional materials with strongly bound excitons, there exists a substantial gap be- 
tween triplet and singlet energies. For example, the lowest triplet (Ti) is about 1 eV lower 
than the respective 5*1 singlet state in conjugated polymers (see Table [11). Our calculations 
show a similar scenario in CNTs owning to quasi-one-dimensional structure of the material. 
In the (7,6) tube there is 0.3-0.4 eV gap between Ti state and optically allowed state. 
This translates into 0.2-0.3 eV gap between Ti state and optically forbidden lowest 'dark' 
Dll singlet state. We expect that the observed triplet-singlet gap in CNTs will reduce with 
increase of the tube diameter, showing behavior similar to the exciton binding energy scal- 
ing. The delocalization pattern of Ti in the (7,6) tube is strikingly similar to that of the PPV 
oligomer (see Fig. [2]). Ti state in both materials is a tightly bound exciton with maximal 
separation between an electron and a hole not exceeding 1 nm. Our results dispute previous 
theoretical study based on the empirical model approximation, which assigned triplet states 
to be within 20 meV of the singlet bands^^. The latter prediction represents a typical solid 
state case and contravenes strong exciton binding energy found in CNTs^i^i^i^ii^. 

Distinct triplet absorption bands provide a major experimental tool for detection and 
understanding of triplet states--. Similar to the singlet manifold with optically allowed 
Si — Sn band, there exists well defined Ti — transition in both polymers and nanotubes. 
Tn states have long been explored in the luminescent polymer's context (e.g.— i^), where 
Ti — Tn splitting is about 1.4-2 eV. Our calculations predict Ti — T„ splitting for (7,6) tube 
to be about 0.5-0.6 eV (see table [T]). T„ excitation is a delocalized exciton, which has a 
structure very similar to that of Sn transition, as evidenced by the transition density plots 
in Fig. O Indeed, with substantial electron-hole distance, the spin direction becomes of lesser 
importance compared to the tightly bound excitons. Consequently, Sn — Tn separation is 
small compared to the — Ti gap in conjugated polymers. This Sn — Tn difference becomes 
even smaller (50 meV) in CNTs as predicted by our calculations. This gives a powerful clue 
where to look for Tn state in CNTs spectroscopies: right below the two-photon allowed Sn 
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excitation. 

In conclusion, we have investigated singlet and triplet excited state manifolds of carbon 
nanotubes using time-dependent density functional theory. A subset of essential singlet 
and triplet states, that dominate photophysical properties of CNTs, is analyzed using cal- 
culated transition densities. These CNT states have very similar properties compared to 
the analogous states previously observed in conjugated polymers. Good agreement of our 
computational results with available experimental data ensures reliable prediction of triplet 
energies in CNTs. Our results place CNTs in the same category of many molecular materials 
such as acenes, luminescent polymers, and MX chains^Si^i^i, where the energy of the first 
triplet state is typically 2/3 that of the first singlet excited state. These evidence significant 
electronic correlation effects in CNTs. The lowest CNT triplet excitations Ti are spatially 
localized with an excitonic size of about 1 nm. We predict strong optical absorption from 
Ti to the higher lying triplets T„ at about 0.5-0.6 eV. T„ excitations are delocalized states, 
which are energetically slightly below the two-photon allowed state Sn in CNTs. These re- 
sults provide specific guidehnes, which make possible experimental detection of triplet state 
in CNTs. The observed subtle interplay between singlet and triplet manifold energetics 
needs to be accounted for when designing specific light-driven or electronic nanotechnologi- 
cal applications based on carbon nanotube materials. We recall that in conjugated polymers 
triplet states are the dominant species formed on charge recombination which yields elec- 
troluminescence. Substantial deviations from the spin statistics (i.e., one singlet exciton is 
formed for every three triplets) favor singlets and higher luminescence yield, and have been a 
subject of intense recent debate^^. Due to spin statics, formation of triplet excitons from the 
electron and hole carriers is possible in transport processes, for example, in the CNTs-based 
chip structures^i^. As evidenced by our results, the CNT triplet states have lower energy 
and even more tightly bound compared to the singlets. Such excitons can strongly affect the 
device performance. Finally, we point out to an interesting possible application of CNT due 
to low-lying triplet states: photoprotection against triplet states and singlet oxygen such as 
carotenoids functions in photosynthesis^. 
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TABLE I: Calculated and experimental excitation energies (eV) of poly-phenylene-vinylene (PPV) 
and (7,6) tube. Experimental values have been reported in ^iS^iSI^ 'ii^ "liiSi. 



Compound 


Transition 


B3LYP 


PBEIPBE 


B3LYP 


PBEIPBE 


Experiment 






ST0-3G 


ST0-3G 


6-3 IG 


6-31G 




PPV 


So — Si 


2.80 


2.99 


2.24 


2.37 


2.48'' 


16 units 


So — Ti 


1.77 


1.64 


1.53 


1.44 


1.3'' 




So — Sn 


3.14 


3.42 


2.59 


2.85 


3.2'' 




So — Tn 


3.05 


3.32 


2.50 


2.75 


3.0" 


(7,6) 


So — Si 


1.39 


1.48 


1.24 


1.29 


l.ll'' 


2 units 


So — Ti 


1.05 


0.95 


0.97 


0.94 






So - Dll 


1.29 


1.37 


1.15 


1.20 






So — Sn 


1.65 


1.82 


1.41 


1.47 


'-l.l'^ 




So — Tn 


1.63 


1.78 


1.39 


1.45 
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FIG. 1: Scheme of the excited state structure of vr-conjugated polymers and semiconductor carbon 
nanotubes. Shown by grey arrows are optically allowed electronic transitions and typical related 
energy numbers corresponding to poly-phenylene-vinylene (PPV) and 1 nm diameter tube. Shaded 
area denote the beginning of a continuum for each excitonic manifold. 
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FIG. 2: Analysis of transition density matrices corresponding to optically active excited states 
of conjugated polymers (left column) and carbon nanotubes (right column) calculated at 
PBE1PBE/6-31G level. These states are schematically shown in Fig. [TJ Each plot depicts prob- 
abilities of an electron moving from one molecular position (horizontal axis) to another (vertical 

axis) upon electronic excitation. The color code is shown in the middle. 
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